Spectroscopic properties of inhomogeneously broadened spin ensembles in a cavity 
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In large ensembles of identical atoms or spins, the interaction with a mode of the electromagnetic 
radiation field concentrates in a single superradiant degree of freedom with a collectively enhanced 
coupling. Given a controllable inhomogeneous broadening, such ensembles may be used for multi- 
mode storage of quantum states of the radiation field with applications in quantum communication 
networks and quantum computers. In this paper we analyze how the width and shape of the inhomo- 
geneous broadening influence the collective enhancement and the dynamics of the cavity-ensemble 
system with focus on the consequences for the ensemble's applicability for quantum information 
processing tasks. 
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I. INTRODUCTION 

In ensembles of a large number of identical atoms, the 
interaction with the electromagnetic radiation field con- 
centrates in a few collective degrees of freedom. Typ- 
ically, one identifies effective oscillator degrees of free- 
dom, which represent the collective atomic population of 
different internal states. The strong effective coupling of 
these oscillators to incident quantum fields makes atomic 
ensembles a prospective component in light-matter inter- 
faces [1] , quantum memories [2] , repeaters for long-range 
quantum communication J3^, and many other applica- 
tions in quantum information technology. 

The ensemble size on the one hand provides a large 
optical depth, and on the other hand it provides phase 
matching conditions to couple strongly to weak quantum 
fields. In a microscopic quantum formulation, an incident 
single photon couples to a single collective excitation, i.e., 
to a superposition state where all atoms have the same 
or similar excitation amplitude. The corresponding Rabi 
frequency is proportional to the square root of the num- 
ber of atoms. While a single collective oscillator degree 
of freedom is coupled to the field, a stored field state can 
be transferred to another spatial mode of collective exci- 
tation by applying a controlled reversible inhomogeneous 
broadening (CRIB) to the atoms. This way, a multimode 
optical interface and memory can be established j4j. Al- 
though optical transitions of rare earth ion dopants in 
crystals are inhomogeneously broadened, narrow spec- 
tral features can still be defined with hole burning tech- 
niques, and both CRIB schemes [5j|6] and schemes based 
on atomic frequency combs [7] have been used to demon- 
strate storage of up to 1060 pulses of light [Hlin]- 

In this paper, we consider an ensemble of N effective 
spin-1/2 particles, coupled to a central quantum oscilla- 
tor (cavity) . A cavity can be used to significantly enhance 
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the atom-light interaction, as demonstrated by experi- 
ments on ultracold rubidium atoms [TU] and ion Coulomb 
crystals }llj in an optical cavity. Furthermore, a single 
quantized field mode can be used as an interface in hybrid 
proposals for quantum computing, e.g., to couple polar 
molecules [121 lU or solid state spin ensembles [T4| to a 
superconducting qubit via a transmission-line resonator. 
Recent experimental breakthroughs have led to the obser- 
vation of strong collective coupling between a supercon- 
ducting transmission-line resonator and large ensembles 
{N > 10^^) of electron spins of chromium ions in ruby [T3] 
or nitrogen- vacancy centers in diamond [16] . while iron 
nuclei embedded in a low-Q planar cavity have been reso- 
nantly excited by synchrotron radiation to a superradiant 
state with a large collective Lamb shift [T7]. The multi- 
mode capacity of ensembles has also been observed with 
nitrogen electron spins in fullerene cages and the electron 
and nuclear spins of phosphorous in silicon \18\ . 

The purpose of this paper is to investigate collective 
enhancement in the presence of inhomogeneous broaden- 
ing, which is not controllable by the experimentalist. We 
study a system with an interaction free Hamiltonian 

Ho ^ojcalac + ^uj-jd-l, (1) 
j 

where lOc is the angular frequency of the central oscillator, 
ujj is the time-independent transition frequency of the 
jth spin, and H — 1. The spins are non- interacting and 
the coupling to the central oscillator is described by a 
Jaynes-Cummings Hamiltonian, 

Hi=Y,9j^ia^ + ii-C; (2) 

j 

with possibly different coupling constants gj . 

This model has fundamental importance and it ap- 
pears with different variations in quantum physics: a 
single spin coupled to a bath of oscillators (the model 
explaining spontaneous emission of light from a single 
atom), a central spin coupled to a spin bath in a spin- 
star configuration [TS] , a single oscillator interacting with 
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a spin ensemble, or a cloud of atoms in an optical or 
microwave cavity. Many of these models have been ad- 
dressed in textbooks [20] with emphasis on the resulting 
dissipative dynamics of the central system. The math- 
ematical difhculties arising when taking the limit of a 
continuously dense ensemble are also well-known [2TH23] . 

Here we revisit the problem with emphasis on the "spin 
bath" degrees of freedom. We show that the density of 
spin states plays an important role in the joint dynam- 
ics of the central oscillator and the collective spin wave 
mode that is directly coupled to it. In certain cases, de- 
spite the inhomogeneity, all the other spin wave modes 
are effectively decoupled. This fact manifests itself in 
reduced oscillator linewidths determined by the homo- 
geneous linewidth of the individual spins [24l [25] . For 
other configurations, however, the oscillator linewidths 
are dominated by the inhomogeneous broadening, and 
the "decoherence" of the superradiant mode is collec- 
tively enhanced. 

The paper is organized as follows. Sec. [ll] derives some 
important properties of the model system, including the 
distribution of eigenenergies, the transmission spectrum 
through the cavity, and a formal solution for the dynam- 
ics of the system in the limit of high polarization. In 
Sec. |III[ we analyze in detail the cases where the effec- 
tive density of spin states is Gaussian and Lorentzian, 
and compare in the two cases the oscillator linewidths, 
as well as the degree of Rabi splitting. In Sec.|IV[ we con- 
sider the regime of strong coupling and study the protec- 
tive effects of the Rabi splitting in a perturbative manner. 
Section \V\ summarizes our results. 



N 
3 = 1 



la. 



(5) 



where = "YjCtjCij, defined by the normalized vector 

'^j — 9j/Cl2k IfffcP)^^^! is the creation operator of a delo- 
calized spin wave mode, the so-called superradiant mode. 
This is the concentrated degree of freedom which is cou- 
pled to the cavity with the collective coupling strength 
^ = (Sfe IfffcP)^^^ that scales as -s/iV. An important 
consequence of the Holstein-Primakoff approximation is 
that the excitations become independent (noninteract- 
ing) quasiparticles. Therefore, the dynamics of a single 
excitation provides the general solution, even if the to- 
tal number of excitations is actually much larger than 
unity (but still much smaller than TV, see [32] for the 
case when the number of excitations in the superradiant 
mode is comparable to the number of spins). 

The system consisting of a central oscillator interact- 
ing with a discrete or continuous bath of oscillators as 
described by Eqs. J4| and ([s]) has been studied in many 
physical contexts |l5H23j . In the following, our aim is to 
gain a detailed understanding of the spectroscopic signa- 
ture of the bath, as observed through the cavity mode. 

in Sec. [11 Aj describe the 
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use the resolvent for- 



In order to do so, we wil . 
eigenbasis of H and in Sec. 
malism 20J to calculate the matrix elements of the time 
evolution operator. We find that the eigenbasis as well as 
the evolution operator is closely linked to the level-shift 
function, and we devote Sec. jll C| to studying the analytic 
properties of this function before applying our analysis to 
describe transmission spectroscopy in Sec. jllDj 



II. DRESSED ENSEMBLE 

At low temperatures (high polarization), a collection 
of two-level systems and a collection of oscillators be- 
have identically. This fact is conveniently described 
in the Holstein-Primakoff approximation by introducing 
the bosonic operators cij for each spin-1/2 particle 



^i = -5+«l«J' ^i = a]^Jl-a]aj^a]. (3) 

The nonlinearity introduced by the square root term in 
Eq. ^ ensures that no two excitations can take place 
at the same spin. If we consider delocalized spin waves 
involving a large number of spins compared to the num- 
ber of excitations, the probability that a given spin is 
excited is inversely proportional to the number of spins 
N. Therefore, as long as only a few delocalized spin ex- 
citations are considered, it is reasonable to neglect the 
square root term in Eq. (|3|. In this regime, the free and 
interaction Hamiltonian Q and ([2| become quadratic in 
the bosonic operators (apart from an omitted c-number) , 

N 

Ha = uJcalac + ^ i^ja]aj, (4) 



A. Exact solution of the eigenvalue problem for 
discrete systems 

In the high polarization (low excitation) limit, the total 
Hamiltonian H = Hq + Hi is quadratic in the creation 
and annihilation operators, H = 'Y^ H^^a^^fl^ with /i, 

= c, 1, . . . , iV. We introduce the "dressed" eigenmodes 

^l = Vqcal + ^riqja], (6) 

3 

that bring the total Hamiltonian into the diagonal form 
H = J2q ^q^q^q- I^or atomic spins in an optical cav- 
ity, the quasiparticles corresponding to the excitations of 
these eigenmodes are called cavity polaritons. We will 
refer to rjqc as the photonic amplitude of the qth polari- 
ton mode, while r]gj (j = 1, . . . , TV) is the mode function 
of the spin wave part. In field theory, rjqc is referred to 
as the wave function renormalization constant [23j . 

In what follows, we will assume that the spins are non- 
degenerate. Otherwise, if to > 2 spins have the same 
transition frequency ujj , we can replace them with a sin- 
gle effective spin a'j 9kO,k/gj and an enhanced cou- 
pling constant g'j = (X)™ ISfcP)^^^- The other to — 1 or- 
thogonal combinations are uncoupled. 
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FIG. 1. Schematic graph visuahzing the solutions of Eq. ([7|. 
The intersections of the rational function K{E) (solid) with 
the y = E — ujc hue (dashed) give the eigenenergies Eq. 



The N + 1 eigenenergies in the single-excitation sub- 
space are the solutions for Eq of the equation 



E„ 



- K{Eq) = 0, 



where K is the level-shift function EO] 



K{z) ^ 



1.9,1 



(7) 



(8) 



which we will discuss in more detail in Sec. |II C] The poles 
of the rational function K{z) are located at the different 
values of ujj . For real w^-s, the solutions of Eq. ([7| can 
be found graphically as the intersections of K{E) with 
the line y = E — ojc (see Fig. [T]) . For most of the eigenval- 
ues the perturbation only causes a small shift: Except for 
two modes, E^ < w,„in and i?+ > Wmax, all the perturbed 
eigenvalues are within the region [wmin, Wmax] ■ In fact, as- 
suming increasingly ordered spin transition frequencies, 
wi < W2 < . . . < LOpf, we see that u!q < Eg < w^+i with 
q = l,...,N-l. 

The photonic amplitude and the mode function of the 
normalized eigenmode q can be obtained, respectively, as 



Vqc = 1 



Vqj = 



9jVqc 
Eg - UJj 



(9) 



As an example, Figs.[2^-c show the eigenenergies and the 
photonic amplitudes of the eigenmodes for an ensemble 
of 25 spins in the regimes of strong, intermediate, and 
weak coupling. 



B. Matrix elements of the time evolution operator 

We will now calculate the transition amplitude for 
an excitation created at t = in the oscillator mode 
v to end up in oscillator at a later time t, that is, 
we calculate the matrix elements of the evolution opera- 
tor U{t) = exp {—iHt) on the single-excitation subspace. 
These transition amplitudes are also called Green's func- 
tions and can be expressed as commutator expectation 
values in the Heisenberg picture. 



(10) 



To account for photons leaking out of the cavity at a 
rate k and spins decaying at rate 7-,-, we start from the 
Heisenberg-Langevin equation of motion for the annihi- 
lation operators in the absence of external driving field, 



/c 



(11) 



(12) 



where fc and fj are Langevin noise operators satis- 
fying {mflit')) = 2n5{t-t') and (mfUt')) - 
jjSjkS{t ~ t'), while the cavity frequency and the spin 
transition frequencies are considered complex with neg- 
ative or zero imaginary parts, Imwc = — k and Imcoj — 
— ^7j. For the Green's functions (10) we obtain a system 
of differential equations. 



(13) 



»^G,,(0 = io,G,,it) + g,G,At), (14) 

with the initial condition G^^{0) — 5^u. Substituting 
into Eq. ( 13 ) the formal solution 



Gjc(t) = Gje(0)e-'"^* - ig, / G,,(r)e-'"^(*-^) dr. 



(15) 



we get a closed integro-differential equation for Gcc{t)- 

d r* 

ij^Gcc{t)^UJcG,c{t)-i j Gcc{T)K{t-T)dT, (16) 

where we have introduced the memory kernel function 
i^(t)^5]|5,fe--^*. (17) 



Equation ( 16 1 can be solved in the Fourier domain by 



extending the integral to a proper convolution. For this, 
we introduce the advanced (-I-) and retarded (— ) versions 
of the Green's functions and the memory kernel function. 



G±,(i) ^ T^e{±t)G^,{t), 
K^{t) = TiQ{±t)K{t), 



(18) 



where Q{t) is the Heaviside step function. With these 
functions, we have a proper convolution, 



zjGtAt)=m+^cGtAt)+ / Gt{T)K^{t-T)dT. 



(19) 



Now taking the Fourier transform of Eq. ( 19 ) we obtain 
wG± (w) = 1 -f WcG± (c^) + G± (w)X±(w), (20) 
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FIG. 2. (Color online) Numerical calculation for N = 2h spins with an inhomogeneous width of Aw, sampled from a Gaussian 
distribution. The individual spins decay at a rate 7^ = O.IAlj and are uniformly coupled to the cavity. The cavity decay rate is 
K — O.OSAoj. (a)-(c) The real part of the exact eigenenergies as function of the cavity detuning. Color red (blue) indicates that 
the corresponding polariton is photon-like (spin-like), (d)-(f) Transmission spectrogram, the color hue (saturation) corresponds 
to the phase (logarithmic magnitude) of the susceptibility (38 1. In the strong coupling regime, Q, = 4Auj, (a) and (d), only 

Aid, (b) and (e), a band 
In the weak 



the two extremal eigenmodes have photonic attribute. In the regime of intermediate coupling, 
of intermediate eigenmodes has significant cavity content, and the cavity can "decay" into any of these modes, 
coupling regime, Q = Aa;/4, (c) and (f), there is no collective effect, and only a few, resonant spins are dressed by the cavity. 



from which the expUcit form of the reduced dynamics of 
the central oscillator follows by simple algebra, 



(21) 



The Fourier transforms ^^^{uj) = J Gjj^(i)e' 
conventionally called forward and backward propaga- 
tors. We note, however, that the backward propagators 
G~j^{uj) are undefined if any of the spin relaxation rates 
7j or the cavity k, is non-zero. 

The remaining propagator matrix elements can all be 
expressed in terms of G'^^{uj): Introducing fcj(i) 



and defining fc- [t) similarly to Eq. (18), we find by 



Eq. ( 15 ) that 



- gl~k^{oj)GtM, 



(22) 
(23) 
(24) 



Furthermore, for the amplitude of transiting from the 
cavity to the superradiant spin wave mode, Gsc{t) = 
{[h{t),al{0)\), we have 

G± (^) = GtM = n-'K^iuj)Gti^), (25) 
G±(^) = l-if±(c.)G±(c.). (26) 

C. Analytic properties of the level-shift function 

We shall see that the complex analytic extension of the 
level shift K^{uj) is of special importance [20]. Here we 
summarize some of its properties. First we note that the 
Fourier-Laplace transforms k^{z) = J k^ (t)e''^* dt are 
defined on complementary halves of the complex plane, 



I undefined, if Im z ^ Im 



if Im z ^ Im Uj , 



(27) 
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FIG. 3. The branch cut of the level-shift function K(z) dis- 
played as the straight line uj — |7hom is^ G K) in the complex 
plane. K(z) has a jump of K~ (z) — K^{z) — iVdz) when 
passing through the cut from above. Depending on whether 
K is smaller than, equal to, or larger than |7hom, the pol es o f 
the cavity-cavity propagator, i.e., the roots E± of Eq. (39 1, 



may lie above the cut, on the cut, or on the second Riemann 
sheet below the cut, respectively. 



and when Im z approaches Im ujj — 



from above or 



below, k^{z) and fc - (z) tend to different distributions. 



lim fcf (w- |7j iiry) =V- 



Rew, 



=p iTrS{uj — Re 



(28) 



where uj is real and V denotes the principal value. 

The whole ensemble may consist of discrete and con- 
tinuous sets of spins, and K^{z) — J2j Idjl'^^j^i^) is de- 
fined for Im^; > sup^ Imwj, while K~{z) is defined for 
Im z < vcdj Im Uj . For simplicity, we assume that the 
decay rate is the same for all spins, 7^ = 7hom- If the 
coupling density profile p{u]) = '^^\gj\'^5{uj — Yleujj) is 
an analytic function of the real variable w, then 



k{z)^ 



duj 



(29) 



is a complex analytic extension of K^{z), that has pole 
singularities in the discrete points ujj and a branch cut 
discontinuity along the line uj — |7hom, such that it ap- 
proaches different values from above and from below. 



lim Kiuj 



j7ho 



± iT]) = K"^ (u 



^7ho 



(30) 



as illustrated in Fig. [3j The real and imaginary parts 
of K^{z) = Ac(z) =F ^^c{z) on the branch cut line z 
w — |7hom follow from Eqs. (28) and (29), 



Ac(w - |7hom) 



|7lK 



2'Kp{u)). 



(31) 
(32) 



Thus the analytic continuation of K'^{z) is given by 
Eq. ( 29 ) above the branch cut, while it explores the sec- 



ond Riemann sheet of K{z) below the cut, 



K+{z) = 



K{z) 

A,(z) - |r,(z) 
k{z) - iT,{z) 



if Im z > — 
if Im z — ~ 
if Im z < — 



2 Thorn 1 
2 Thorn 1 
2 Thorn • 



(33) 



We can revert Eq. ( 32 1 and extract information about 
the coupling density from the values of the level-shift 
function j^+(aj) which is evaluated on the real axis, 



1_ 



-linK+[uj - |7hom) 

1 5Re^+(a;) 



lmK+{u}) 



2tt 



-7h< 



(34) 



where we truncated the Taylor series of Im_ftr+(z) at 
first order in 7hom and obtained the derivative from the 
Cauchy-Riemann equations. 

Finally, we remark that we included 7hom in Eq. ( 29 ) 
merely in order to separate the effects of homogeneous 
and inhomogeneous broadening. Alternatively, we could 
have assumed a coupling density profile that already in- 
cludes homogeneous broadening via the convolution 



7hom/27r 



p{uj')duj'. (35) 



The branch cut of the new level-shift function K'{z) = 
J jJ (lS)I(z — uj) duj is then on the real axis. However, the 
analytic continuations of K'^{z) are the same as those of 
K^{z)^ and the two models lead to the same dynamics. 



D. Linear response and spectroscopic features 

An experimental tool for studying the eigenenergies of 
the coupled cavity-ensemble system is transmission spec- 
troscopy. In this process, the (lossy) cavity is driven by 
an external classical field of frequency uq via the inter- 
action Hamiltonian V{t) = f e~*"°*ajj -|- H.c. The trans- 
mitted field — which is proportional to the steady state 
of the field inside the cavity — reflects the linear response 
of the system to the driving field and shows a resonance 
when the driving frequency is near an eigenfrequency. Al- 
ternatively, one may consider, e.g., directly driving the 
spins via an external field that couples to the total spin. 



Vit) = Et 



-icjQt 



H. 



Let us consider a general external classical drive de- 



— scribed by V{t) 



"*"o*a|, -f H.c. We assume that 



the driving field is weak, so that the total number of 
excitations is always small and the Holstein-Primakoff 
approximation remains valid. Then the Heisenberg- 
Langevin equation of motion for the annihilation oper- 
ators can be obtained from Eqs. ( 11 ) and ( 12 ) by adding 

-icjjQt 



i£r 



-iujQt 



and 



-iSjC 



respectively, to the right- 



hand sides. The formal solution of this set of inhomo- 
geneous differential equations can be written in terms of 
the Green's functions as 

a^.{t) = Y,G^,{t~to)aAto) + J2 f Gf^^t-r) 

X [-zf,e-'"»^ + /,(T)] dr. (36) 
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The first term vanishes when taking the limit to ~^ — oo, 
for all eigenvalues of if^^ are assumed to have negative 
imaginary part (even if infinitesimal). Since (/^^(t)) — 0, 
the expectation value of the oscillator operators in the 
steady state becomes 



(37) 



the limiting cases of infinitely narrow and infinitely wide 
coupling density profiles, before studying the qualitative 
differences between the intermediate cases for Lorentzian 
and Gaussian coupling density profiles. 



A. Oscillation and decay 



where 

/>oo 

X^A^) = -I / G^.(r)e*"- dT - G+,(u;) (38) 

JO 

is the susceptibility (impedance) of the system, which has 
been calculated in Eqs. (21)-(26). 



In the transmission spectroscopic setup, the cavity field 
is {a^cit)) — £e~'^'^°*Xcc{^o), so the transmissivity of the 
system is proportional to |G+ (wo)^, and the phase shift 
is 0(wo) = argG'+(a;o) [13 HZ]. The resonance peaks 
in the transmission spectrum are characterized by the 
complex poles of G+(z), that is, the solutions for z of 
the implicit equation 



K+{z) = 



(39) 



For discrete systems, the analytic continuation of K^{z) 
has pole singularities in the points Uj in the lower half of 
the complex plane, and the solutions of Eq. pO] ) coincide 
with the eigenenergies Eq obtained form Eq. ([t]). For 
continuous systems described in Sec. II C the roots of 
Eq. (39) may be located above the branch cut of K{z), 
on the cut line, or below it depending on whether k or 
|7hom is larger (see Fig.|3|. 

Finally, we mention that there is a direct way to 
extract information about the coupling density profile 
from a measured transmission spectrum. First the level- 
shift function is to be calculated by inverting Eq. (21 1: 
K~^{uj) — u! — ojc — x7ci'^)- Alternatively, if only the 
transmissivity is available experimentally but not the 
phase of the susceptibility, K~^{uj) may be obtained by 
fitting a Lorentzian to the transmissivity |xcc('^,Wc)P as 
function of the cavity frequency lUc (with Imwc — --ik). 
Once the values of the level-shift function K^{uj) are 
known for real frequencies, Eq. (34) can be used to derive 
the coupling density. 



III. EXAMPLES OF SPECTROSCOPIC 
SIGNATURES 

With the results of the previous section, we are 
equipped to predict the spectroscopic signature for a 
given system as observed by transmission spectroscopy. 
While the output signal could have been obtained more 
directly via classical input-output theory [25], we shall 
see that the framework developed in Sec. [H] provides ad- 
ditional insights into the physics by explicitly consider- 
ing the state of the spin bath. We will first consider 



In certain situations, the detailed structure of the cou- 
pling density profile is not important. On time scales 
much smaller than the inverse inhomogeneous width of 
the ensemble, that is, when the inhomogeneous broaden- 
ing cannot be resolved, the ensemble behaves as a sin- 
gle oscillator with a collectively enhanced coupling con- 
stant. In this case, the memory kernel function (17) can 
be approximated by K{t) k, fi^e"*'^"*, or equivalently, 
k^iuj) = Ac(w) T |rc(a;) with Ac(w) « V^'^/{uj - w^) 
and Tc(io) « 2tt^I'^6{lu — uja), where the ensemble's mean 
transition frequency uja is assumed to be real. The inverse 
Fourier-Laplace transform of Eqs. (21 1 and (25) show 
that the cavity and the superradiant mode undergo Rabi 
oscillation at the Rabi frequency fij^ = yjW' + (5^/4, 

Gc.c{t) = e"*("°+''/2'* ( cos ^Rt - i cos 61 sin ^^1) , (40) 



Gsc{t) 



- _,V-*("o+'5/2)t 



sinflsinriflt, 



(41) 



where the detuning of the cavity, 5 — lJc — uJai and the 
mixing angle tanfl = 2^/5 are assumed to be real. 

In the opposite limit, when the inhomogeneous width 
is much larger than the collective coupling VL, and the 
spin transition frequencies are sufficiently dense to form 
a continuum, we are in the Weisskopf-Wigner regime. 
In analogy with the spontaneous emission of an excited 
atom, the cavity excitation effectively decays into the 
continuum of spin waves in this regime j28j . The key 
assumption in the Weisskopf-Wigner approximation is 
that the kernel function K{t) decays fast compared to 
the slowly varying envelope Gcc{t)e^^''^ , and thus the fre- 
quency dependence of K^{oj) can be neglected. The in- 
verse Fourier-Laplace transform of Eq. (21) then gives 
an exponential decay, 



Gccit) 



_ --iK+A,(we)-^r,(t^,)]t 



(42) 



Here we see that Ac{ujc) is the shift of the cavity fre- 
quency and ^Tc{ci>c) is an additional cavity decay rate 
due to the coupling to the dense spin reservoir. 

Finally, we consider the general case when an initially 
photonic excitation decays into the continuum of spin 
wave modes, and we specify the asymptotic energy dis- 
tribution of the excited spin. If Gccit) converges to zero, 
it is either because the initially photonic excitation leaks 
out of the cavity, or because the photon is irreversibly 
converted into a spin-like excitation. From Eq. (15) we 
have 



lim G 



(43) 
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so the asymptotic probability that the jth spin is excited 
equals limt_j.oo \Gjc{t)\'^ = assuming that 

jj =0. Then the asymptotic energy distribution of the 
excited spin, i.e., the probability that the transition fre- 
quency of the excited spin is w, reads 



j 



(44) 



In particular, it is proportional to a Lorentzian in the 
Weisskopf-Wigner regime, in analogy to the energy dis- 
tribution of a photon spontaneously emitted by an atom. 



B. Lorentzian coupling density profile 

When the cavity is coupled to a single atomic oscilla- 
tor (e.g., a degenerate ensemble of homogeneously broad 
but identical atoms) , we recover the well-known problem 
of a driven two-level atom or that of two coupled oscil- 
lators. In this simple case, the results of Sec. |II A| can 
be directly applied. With an atomic transition frequency 
Wq, natural atomic linewidth 7, cavity detuning from the 
atomic transition (5, and cavity linewidth 2k, we take the 
complex frequencies lji ^ uja — |7 and = + 5 — ik. 
The eigenenergies obtained from Eq. are 



I7 + f^T^ COS^ ± 



where 



(45) 



(46) 



is the (complex) Rabi frequency, and for the (complex) 
mixing angle we have 



cost 



5 — i{K 2 



sinfl 



(47) 



In analogy with the dressed states of a driven two-level 



atom, the eigenmodes are $j_ 
photonic amplitudes 



r]±cal. + 'r]±sW with the 



V+c 



cos(6'/2) 
■\/cosh(Im6') ' 



V-c 



sm{0/2) 
■\/cosh(Im(?) ' 



(48) 



and the superradiant amplitudes ri±s — ±?7=pc 



The memory kernel function defined in Eq. (17) is 



K{t) = rPe "^°* ^'/2 for such a single collective spin ex- 



citation, which yields K^{u!) — Q,'^/{uj — i 



h)- The 



same memory kernel function is obtained for a large, in- 
homogeneously broadened ensemble with a Lorentzian 
coupling density profile. This is the case, for example, 
when each spin in the ensemble is uniformly coupled to 
the cavity with the same coupling constant g = f2/\/]V, 



and the spin transition frequencies are distributed ac- 
cording to 



2?H = 



1 



iV7 



27T (c^-C^,)2+^2/4' 



(49) 



where now 7 is the ensemble's inhomogeneous width (full- 
width at half-maximum) and the individual spins do not 
decay, and the coupling density profile is p(w) = g^V^co). 
All the dynamical properties of both the cavity and the 
superradiant spin wave mode are determined by K^{uj) 
alone and, thus, they are the same for a single (homo- 
geneously broad) atomic oscillator and for an inhomoge- 
neous ensemble with Lorentzian coupling density. There- 
fore, it cannot be determined from the linear response of 
the cavity whether the broadening is of homogeneous or 
inhomogeneous origin. 



The eigenenergies (45) are the two complex poles of 
the forward propagator (21), 



[uj - E+){uj - E_y 



(50) 



corresponding to two peaks in the transmission spectrum: 
The position and the width of the peaks are determined 
by the real and imaginary parts of E± , respectively. The 
system undergoes damped Rabi oscillations, as can be 
seen, e.g., by taking the inverse Fourier transform of 
Eqs. (l50| and (pSl), 



Gcc(t) = ^-i{'^a+S/2)t-(-,/2+n)t/2 

( ^ cos Re Q 

X cos iZi^i — i- 



sinrfflt , (51) 



cosh Im Q 

X (sinRc(?sinilflt — tanhIm6'cosr2flt). (52) 

As an illustrative example, let us consider the resonant 
case 6 = 0. Depending on the value of (7 — 2k) /(4ri), we 
see a transition from (i) damped Rabi oscillations at a 
reduced effective Rabi frequency to (ii) an overdamped 
situation (see also Figs.|4^-b). 

(i) Regime of oscillations. The Rabi frequency (46 1 is 
real when 7 — 2k < AH., giving rise to a Rabi splitting of 



y^rj^ _ _ k)2/4 around the atomic transition 



frequency uja (Fig. |4p,) , and the width of the peaks in the 
transmission spectrum is — 2Im£'-|- = ^7 -I- k (Fig. [ija). 
Here we see two important facts about a Lorentzian 
coupling density profile: the resonant Rabi splitting is 
smaller than the actual collective coupling strength il, 
and the linewidth of the peaks is determined by the in- 
homogeneous width of the ensemble, independent of the 
coupling strength. As we shall see, the cavity does not 
provide protection for the superradiant mode in this case. 

(ii) Overdamped regime. The Rabi frequency (46 1 is 
pure imaginary for 7 — 2k > 4f2. There is no Rabi split- 
ting at all, for we have two overlapping peaks at uJa of 
width -2Im E± = {j + 2k) /2 =F ^(7 - 2k)2/4 - 41^2. In 



8 



Gaussian 




0.0 0.2 0.4 0.6 0.8 

^/ 'I'fwhm 



0.2 0.4 0.6 0.8 1.0 

^/Tfwhm 



FIG. 4. The real (a, c) and imaginary (b, d) parts of the 
poles of the susceptibility (21 1 as function of the collective 



coupling strength for a Lorentzian (a, b) and a Gaussian (c, d) 
coupling density profile. The quantities are directly related 
to the position and width of the peaks in the transmission 
spectrum. The cavity is tuned at the center of the ensem- 
ble (ujc = ^a), and K — 7hom ~ was assumed. The full- 
width at half-maximum of the coupling density is 7fwhm ~ 7 
for the Lorentzian and 7fwhm ~ VS In 2 ct for the Gaussian. 
The dashed lines (a, c) show the lowest order solutions of 
the asymptotic expansion given by Eq. (77 1, and correspond 
to the strong coupling approximation described in Sec. [Ill K\ 
The dotted curves (c, d) show the asymptotic solutions in the 
next orders as given by Eqs. (62 1 and (63 1. The dash-dotted 



curve (c) shows Eq. ( 58 1 corresponding to the Weisskopf- 
Wigner approximation in the weak coupling regime. 



the Weisskopf-Wigner regime (7 ^ ft, 2k), the dynamics 
reduces to the exponential decay given by Eq. ( 42 ) , and 



the peaks in the spectrum become clearly distinguishable: 
On top of a wide background peak with 7^ « 7 corre- 
sponding to the weakly perturbed superradiant mode, 
there is superimposed a narrow peak of the cavity, whose 
width 7_ sa 2k -t- 4ri^/7 equals the overall decay rate of 
the cavity excitation as given in Eq. ( 42 ) . 



Finally we note that the oscillating-decaying dynam- 
ics of our cavity photon interacting with an inhomoge- 
neously broadened, highly polarized spin ensemble shows 
some similarity with the dynamics of a single atom in a 
damped, zero temperature cavity |29j . In the latter case, 
the atom interacts with the lossy cavity, and the cav- 
ity itself is coupled to a collection of external radiation 
modes. In our case, the cavity interacts with the superra- 
diant spin wave mode, and the superradiant mode itself 
is coupled to a collection of subradiant spin wave modes. 



C. Gaussian coupling density profile 



We consider now a Gaussian coupling density profile. 



p{uj) 



1 



2tt cr 



exp 



{UJ - UJaY' 
2cr2 



(53) 



and show that the spectrum of an ensemble with such 
an inhomogeneity is qualitatively different from what we 
learned from the Lorentzian profile in the previous sec- 



tion. From the integral in Eq. (291 we obtain 



K+{z) 



2 a 



erfi(0 



(54) 



where we have introduced the dimensionless parameter 
I + 57hom) /(VScr). Equation (ImI) is com- 



posed of complex analytic functions with no branch cut 
singularities, and Imiir+(a;— |7hom) = —Trp{uj). There- 
fore, Eq. (54) is indeed the analytic continuation of 
K^{z) on the entire complex plane. To determine the 
position and width of the peaks in the transmission spec- 
trum, we solve Eq. (39) numerically (Figs.|4j:-d). Here we 
observe two regimes: an overdamped one without Rabi 
splitting (f2/7FWHM ^ 0.23) and another one with two 
Rabi-split peaks {^/^fwkm ^ 0.23). The most impor- 
tant difference from the case with Lorentzian coupling 
density is that the width of the split peaks decreases with 
the coupling strength, and for strong coupling it is domi- 
nated by the ensemble's homogeneous width or the cavity 
K, as has been pointed out in Refs. [241 125] . 

To estimate the position and width of the resonances, 
let us first consider the overdamped regime, ft a. The 
Taylor series expansion of Eq. ( 54 ) around ^ = reads 



V2o 



(-2^2)" 
(2rH-l)!! 



i-eY' 



(55) 



A good approximation for the first peak, which corre- 
sponds to the dressed cavity mode, can be readily ob- 
tained by keeping only the first term in both sums, 

Re{E+ -uja) = 6/{l-n^/(7^), (56) 
2k - 7hom + V^na 



^2ImS4 



2k 



(57) 



For 6 



K = 7hom 0, (dash-dotted curve in Fig. [4]: 

-2Im£;+ = V2^n^/<7 + 0{n'^) (58) 



has quadratic dependence on the collective coupling Q, 
which confirms the Weisskopf-Wigner approximation. 
Keeping the term in Eq. (55), we find a second root. 



Re {E^ - uja, 
-21vaE- — 27hom ~ 2k 



8a 



2^ Vf^' 



- 1 



(59) 



(60) 
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It is worth mentioning that we have not only two, but in- 
finitely many roots with increasingly large negative imag- 
inary part (not shown in Fig. [4]). This is a mathematical 
consequence of the dynamics of the subradiant modes and 
causes the cavity not to follow a simple exponential de- 
cay |30j . In the limit of weak coupling, however, the first 
root with finite imaginary part yields exponential decay 
in accord with the Weisskopf-Wigner approximation. 

We consider now the strong coupling regime, f2 3> c 
The asymptotic expansion of Eq. ( 54 ) reads 



V2o 



E 



(2n- 1)!! 



2^n 



(61) 



Although the asymptotic power series does not converge, 
for any given ^, we still get a good approximation if we 
truncate the series at n < In fact, the strong cou- 
pling approximation in Sec. |III A| corresponds to keeping 
only the n — term of the series, and it yields eigenener- 
gies similar to Eq. (45) with 7 — 7hom, which are shown 
as dashed lines in Fig, 
(5 = and k = 7hom = 



4p. Focusing on the resonant case 
0, we first neglect the imaginary 
part in Eq. (61) and truncate the series after the n = 1 
term. This gives the position of the peaks in the next 
order (dotted curves in Fig.[4|;), 



Re{E. 



(62) 



or equivalently, = ±^/{n'^ + a'^)/2a^. To obtain the 
imaginary part of the eigenenergies, we look for the solu- 



e. We write ^j."'^' in place of ^ 



tions in the form ^ = ' 

in the exponent in Eq. (61), and we use the facts that 
and 1/C^ 



With these approximations we have for the imaginary 
parts (dotted curve in Fig. |4jj 



2a 



(63) 



IV. PROTECTIVE ENERGY GAP IN THE 
STRONG COUPLING REGIME 

Ideally, when a spectrally narrow ensemble is strongly 
coupled to a cavity, excitations undergo coherent Rabi 
oscillation between the cavity and the superradiant spin 
wave mode at a collectively enhanced Rabi frequency Q, 
without involving the subradiant modes. In reality, how- 
ever, the inhomogeneity in the spin transition frequencies 
will gradually mix in the subradiant modes, thus result- 
ing in decoherence of the cavity-superradiant subspace. 
The time scale for the latter process can vary, depend- 
ing not only on the ensemble's inhomogeneous width, 
but also on the structure of the inhomogeneity. For a 
Lorentzian coupling density profile, e.g., it is always dom- 
inated by the inverse of the ensemble's inhomogeneous 
width, as we have seen in Sec. |III B[ For other distribu- 
tions, like the Gaussian, this "decoherence" time can be 



significantly longer — an effect conventionally explained 
by a gapping mechanism [301 131j . In this Section, we in- 
vestigate how the strong coupling to the cavity prevents 
the mixing of the superradiant and subradiant modes, 
thus leading to narrowing of the line width of the su- 
perradiant spin wave mode. For the cavity to have an 
effect, we will assume that the cavity is not too far from 
resonance {S <C fl^/Auj). 



A. Appearance of the energy gap 



In the absence of inhomogeneity, as shown in Sec. |IIIB| 
the "+" and "— " eigenmodes are combinations of the 
superradiant spin wave mode and the cavity mode, 



= cos I dc -I- sin | b, 
= sin I flc - cos I b, 

p{0) 



(64) 
(65) 



with the eigenenergies i?^ — {uJ + Wc)/2 ± flR. In 
the presence of inhomogeneity, however, these polariton 
modes are only approximate eigenmodes. Their equa- 
tions of motion read 

$f = -j£;f $f - tAiu sin I c, (66) 



dt 



= -iE'^'><^>'^' + i Aw cos I c, 
dt ^ 



r(0)A(0) 



(67) 



where we omitted the Langevin noise terms, c = 
Auj~^ ^j{ujj —uj)a*dj is an annihilation operator cor- 
responding to a subradiant mode orthogonal to both dc 
and 6, and the ensemble's inhomogeneous width Aw is 
defined as the variance of the spin transition frequencies. 



Aw^ = E 



E-.Ki' = / 



j7hom) du, 



(68) 
(69) 



where we assumed Imw, = Imw 



uj — iiiiv^ — — ^7hom for all spins. 
For a cavity not too far from resonance {6 <C ri^/Aw), 
the limit of strong coupling is identified by the condi- 
tion n ^ Aw. If the tail of the coupling density profile 
falls off sufficiently fast so that the Rabi-split eigenener- 
gies E^^^ lie far from all the spin transition frequencies, 
then the operators <l>j^''' rotate fast with respect to c, and 



the contribution of c in Eqs. (66) and (67) can be ne- 



glected (rotating wave approximation). In other words, 
the spin dephasing processes induced by the inhomoge- 
neous broadening have to bridge the energy gap between 
the dressed modes and the subradiant modes, and if this 
energy gap is large enough, the dressed modes are effi- 
ciently protected from decoherence (see Fig. [5]). 

Finally, we note that the gapping mechanism may 
be observed even when the cavity is out of resonance 
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FIG. 5. (Color online) Level diagram showing the quasi-closed 
subspace consisting of the dressed polariton modes (superpo- 
sitions of the cavity and the superradiant spin wave mode) 
coupled to a chain of subradiant spin wave modes. In the 
strong coupling regime Qr ^ Auj, the coupling between the 
dressed modes 4>± and the subradiant mode c is off-resonant, 
therefore, the inhomogeneity does not induce decoherence of 
the dressed modes: the dressed levels are not broadened. 



(Aw < < ^ < n'^/Aut). Provided that the tail of the 
coupling density profile falls off sufficiently fast, the cav- 
ity can be adiabatically eliminated in the Born-Markov 
approximation. In this case, the cavity only has a dis- 
persive effect, and its presence results in an energy shift 
of the superradiant mode equivalent to the ac Stark shift 
in optics. Namely, Hi can be replaced by 



(70) 



independent of the state of the cavity. The superradiant 
mode, = ^jOijCL^p is an eigenmode of -ffgap, but due 

to the dephasing effects of the inhomogeneity in Hq, it 
gets mixed with the subradiant modes, as shown by the 
Heisenberg equation of motion 



dt 



6^ =^[i^0+i^gap,&1 -*(i^+Agap)St+zAL^ct. (71) 



For a large gap (Agap ^ Aw) and provided that the tail 
of the coupling density profile falls off sufficiently fast, 
the superradiant mode is energetically separated from the 
subradiant modes. Therefore, the inhomogeneous broad- 
ening cannot induce real transitions involving the super- 
radiant mode. This energy gap may efficiently protect 
the quantum information stored in the superradiant spin 
wave mode from dephasing effects of the inhomogeneous 
broadening or spin diffusion ^31) . 



B. Corrections to the eigenenergies 

In this section, we estimate how fast the tail of the cou- 
pling density distribution should fall off in order for the 
cavity-superradiant subspace to become protected and 
for the line narrowing to manifest. In the strong cou- 
pling limit, the extremal eigenenergies E± lie far away 
from the typical spin transitions frequencies, and in this 
region (|z — a;| ^ Aw) we make an asymptotic power 



series expansion of K^{z) in the form 



fc=0 



o 



1 



[Z - w)"+2 



(72) 



For the Lorentzian (49), An = (— 57)", while A2n = 



{2n - l)!!cr2» and A2n+i = for the Gaussian ([53]). In 
general, inserting the identity 



V- 



1 



00 



Rew)" 



^ (w - Rew)"+i 

n=0 ^ ' 



(73) 



into (311 yields the asymptotic power series expansion 

00 ^ 

Ac(w - f 7hom) ^ ^ J— 



^ (w - Rca;)"+i 



(74) 



as long as the statistical moments of the coupling density, 
1 



[uj — Rew)"/9(w) dw. 



(75) 



are well-defined. In what follows, we will assume that 
the coupling density also has an asymptotic power series 
expansion. 



11 Br, 



^M-E(,„Rew)"+i- 

Tl — 1 ^ ' 



(76) 



Then ImA^ = — Tri?^. However, if the n-th moment (75) 
exists, we must have = 0, and A^ — Mk for all k < n. 



Given the power series expansion ( 72 ) , we now look for 



the resonance peaks in an iterative way. Since Mq = 1, 
we start with _ft'+'-'^^(z) k. Vi^ /[z — uj) in the zeroth order. 
Equation ( 39 ) then yields the same roots as those given 
by Eqs. (|45P(|47l) with 7 = 7hom, namely 



E^^^ - Uj = ±n{cot 



(77) 



In the first order, we look for the solution in the form 

z = E^^'' + e and write 



1 



A, 



E(°) - UJ Ei°) - UJ 
Then we obtain the corrections to S^^' and E^"^ , 

.(1) 



= Ai sin" 



e_ = Ai cos 



2 e 
2 • 



(78) 



(79) 



In the case of a Lorentzian coupling density, Ai = — §7 
immediately leads to wide resonance peaks dominated by 
the ensemble's inhomogeneous width 7. In general, if the 
tail of the coupling density profile falls off as (w — w)~^, 
then the imaginary part of the correction e'^-* reads 



-2Ime 



(1) 
± 



27rBi 



sm 
cos^ 



(80) 
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and it is not reduced by a strong coupling constant fi. If, 
however, the first moment of p{uj) exists (and U is such 
that Ml = 0) and p{uj) has the asymptotic power series 
expansion (76), then Bi must be zero, and so Ai — 0. In 



this case, we have to proceed further in the expansion. 

In the second order, we assume that Ai = 0, so we 
look for the solution in the form z — E^'^^ — e and write 



E(0} - , 



■Ao 



(81) 



To simplify the roots of the resulting quadratic equation, 
we assume that \A2\ <C and make a series expansion 
with respect to ^2/^^- Then we obtain for the corrections 



.(2) 



A| fsin^ltan^, 2/^4, 
n \cos2|cot I / ^ 2/ 



and the resonance peaks are at 



E_ 



(2) _ UJ+UJc 



E. 



(2) 



2 

To + LOr 



(2sin^f)|' 
(2cos^f)| 



(82) 

(83) 
(84) 



We see that the contribution from A2 is suppressed for 
large fi. If the second moment of p{uj) exists, then A2 = 
Aw^ is pure real, and corrections to the resonance widths 
may only come from a complex mixing angle 9. As an ex- 
ample, let us consider the resonant case (5 = 0. The Rabi 
frequency is reduced, flR, = [il^ — (|7hom " k)^/4]"^^^, 
and the mixing angle is « 7r/2 — i(^7hom — K)/2f2fl. 
The real and imaginary parts of the two poles then read 



Re (^i'^ 



1 Aw^ 
2 



^2 



-2Im£;. 



(2) 



llho 



+ K + (^Thora - k) 



Aw^ 



(85) 
(86) 



A small inhomogeneous broadening enhances the Rabi 
splitting, while the width of the peaks is predominantly 
determined by the homogeneous width, and not the inho- 
mogeneous one. We mention that in contrast to Eq. (86 1, 
Eq. (63) was obtained for k = 7hom = and, therefore. 



we had to consider higher orders of p{i-o) to estimate the 
peak widths. 



C. Losses due to inhomogeneous broadening 

To estimate the magnitude of the leakage from the two- 
dimensional subspace of the dressed states due to the 
inhomogeneous broadening, consider an initial excitation 
in the mode ( [64| ) . The exact time evolution of this mode 
in the Heisenberg picture reads 



<l>f(t) 



9' 



(87) 



with = COS ^r]qc+am ^rjqs and q = +, ~, 1, . . . , K — 1. 
We expect that the transition amplitude between the 
original and the evolved modes, given by 



G++(i)^([$f (t),<l 



rotates with an amplitude that is only slightly decreased 
compared to the homogeneous case. Indeed, we can place 
a lower bound on |G++(t)| using the triangle inequality, 



\G++{t)\ > 



^|0,P = 2|<^+P-1. 



Using the method of Sec. IV B we obtain 



V+c 



cos : 



sm 



l-sin^|(2 + sec^|) 



202 



l-htan2|cos6'(H-cos2|) 



Aw' 
202 



(89) 

(90) 
(91) 



and the inequality (89) reads, to second order in Aco/Q, 

.2 



\G++{t)\ > 2|</)+|2 - 1 w 1 - 2Re (sin^ f tan^ f ) 



(92) 



which yields |G++(i)| > 1 — for the resonant 

case ^ = 0, and |G++(t)| > 1 - 4AwVA^ap for the off- 
resonant gapping regime {S <C il'^ /Auj, i.e., Aui ^ Agap). 

The leakage (1— |G++(i)|) is thus bounded from above 
by a quantity proportional to Aa;^/02. A similar bound 

can be derived for a excitation. The bound given by 



Eq. ( 89 ) is exact and valid for arbitrary coupling density. 



We emphasize, however, that these bounds are obtained 
by considering only the inhomogeneous distribution of 
the spin transition frequencies. Other relaxation mech- 
anisms of the individual spins, such as spontaneous de- 
cay, dephasing, or spin diffusion, may still be significant 
sources of losses. 



V. CONCLUSION 

We have studied highly polarized, inhomogeneously 
broad spin ensembles interacting with a single mode of 
a cavity. Using the resolvent formalism, we have shown 
how the transmission spectrum of the cavity depends on 
the coupling density profile and how the coupling density 
can be obtained from the transmission spectrum. 

The strong superradiant coupling provides an energy 
gap for both the cavity and the superradiant spin wave 
modes from the subradiant modes, which may, in certain 
cases, efficiently protect the superradiant polariton and 
decrease its linewidth. We have investigated the crite- 
rion for the appearance of this gapping mechanism in the 
presence of inhomogeneity, and provided corrections to 
the conventionally used picture of a driven two-level sys- 
tem. We have also considered two specific inhomogeneous 
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coupling density profiles: Lorentzian and Gaussian. We 
have shown that the gapping mechanism does not work 
for the former, and the polariton linewidth does not de- 
crease with the collective coupling strength. For the lat- 
ter, however, the gapping mechanism can efficiently re- 



duce the polariton linewidth to a limit depending only 
on the cavity linewidth and the ensemble's homogeneous 
linewidth. In general, we have found that the ensemble's 
coupling density p{uj) should fall off as or faster in 
order for the gapping mechanism to manifest itself and 
line narrowing to take place. 
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